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Abstract 

Conservation equations for the mass, linear momentum and energy 
densities of solitons propagating in finite, infinite and periodic, non- 
linear, planar waveguides and governed by the nonlinear Schrodinger 
equation are derived. These conservation equations are used to deter- 
mine classical force densities which are compared with those derived 
by drawing a quantum mechanics analogy between the propagation of 
solitons and the motion of a quantum particle in a nonlinear potential 
well. 

KEY TERMS: Spatial solitons, nonlinear Schrodinger equation, infinite 
planar waveguides, classical forces 
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1 Introduction 



Soliton propagation has received a great deal of attention in recent years 
because of its ubiquity in many branches of physics, e.g., electrodynamics, 
fluid dynamics, plasma physics, etc., and because of its possible applications 
to distorsionless signal transmission in ultra-high speed and long-distance 
transoceanic telecommunications by optical fibers, soliton lasers, all-optical 
soliton switches, nonlinear planar waveguides, nonlinear transmission lines, 
etc. [1]. For example, the envelope of the electric field in a weakly nonlinear, 
planar waveguide is governed by the nonlinear Schrodinger (NLS) equation 
which can be used to determine the forces acting on solitons [2,3]. 

Based on a quantum mechanics analogy between the NLS equation and 
the motion of a quantum particle in a potential well, Bian and Chan [2] 
and Ramos and Villatoro [3] determined the nonlinear and the diffraction 
force densities on solitons propagating in infinite and finite, respectively, 
weakly nonlinear, planar waveguides. In particular, Ramos and Villatoro [3] 
determined numerically the quantum force densities on solitons in finite, 
weakly nonlinear, planar waveguides subjected to periodic and homogeneous 
Dirichlet, Neumann and Robin boundary conditions at both boundaries of 
the waveguide. 

In this paper, conservation equations for the mass, linear momentum and 
energy densities of solitons governed by the nth-dimensional NLS equation 
are derived. In particular, the linear momentum equation is used to obtain 
the classical force densities which are compared with the quantum force den- 
sities derived in Reference 3. Furthermore, these conservation equations are 
employed to assess the changes in linear momentum that a soliton under- 
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goes in a finite, weakly nonlinear, planar waveguide as it collides with the 
waveguide's boundaries for periodic and homogeneous Dirichlet, Neumann 
and Robin boundary conditions. 

2 Conservation Equations and Classical Force 
Densities on Solitons 

As shown in References 1 and 2, the one-dimensional NLS equation in di- 
mensionless Cartesian coordinates can be written as 

iu t = —u xx — \u\ 2 u, x G V, t > (1) 

where V is the spatial domain, the subscripts t and x denote partial dif- 
ferentiation with respect to time and the spatial coordinate, respectively, 
and u denotes the slowly varying amplitude of the electric field in a weakly, 
nonlinear slab waveguide. 

Equation ([I]) can be written in dimensional form by transforming the 
dependent and independent variables as indicated in References 2 and 3. 

In a recent work [3], the nondimensional NLS equation was used to draw 
an analogy between the propagation of a soliton in a weakly nonlinear, pla- 
nar waveguide and the motion of a quantum particle in a nonlinear potential 
well. This analogy allowed to determine both the quantum energy density 
of and the force densities on solitons propagating in finite, nonlinear, planar 
waveguides. The main results obtained in [3] are repeated here for complete- 
ness and comparison with the classical mechanics ones derived in this paper, 
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and can be written as 



e q (x, t) = Tie {iu*u t } = —Tie {u*u xx + |w| 4 } = e k (x, t) + e v (x, t), (2) 

where e denotes energy density, the subscripts q, k and v stand for quantum, 
kinetic and potential, respectively, Tie denotes real part, and the asterisk 
denotes complex conjugation. 

From the quantum energy density, one can easily obtain [3] 

0\u\^ 

f n (x,t) = -—e v (x,t) = -jfe-, fd(x,t) = -—e k (x,t) = —Tle{u*u xx }, 

(3) 

where / denotes a quantum force density, and the subscripts n and d stand 
for nonlinearity and diffraction, respectively. 

Classical conservation equations and classical force densities on solitons 
may be obtained from the one-dimensional NLS equation as follows. The 
mass (to), linear momentum (M) and energy (E) densities of the soliton can 
be defined as 

to = u*u, M = i (u* x u — u* u x ) , E = \u x \ 2 — -\u\ 4 . (4) 

Using the NLS equation, i.e., Eq. (1), and the above densities, the following 
conservation equations can be easily obtained 

dm DM n 

DM d ( \2_ d2m _ m ^\ 9 ( m 2 AE _ d2m \ 
dt dx \ dx 2 J dx \ dx 2 J ' 

dE Of ,du du*\ .d „ d 
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It must be noted that the classical energy density can be written as the 
sum of potential and kinetic energies using action-angle variables [4]. 

For initial- value problems, i.e., — oo < x < oo, t > and \u\ — > as 
| a; | — > oo, the integrals of the left-hand sides of the above three conservation 
equations over the whole spatial domain are zero. Therefore, the total mass 
or number of particles, the total linear momentum and the total energy are 
constant and coincide with the first, second and third invariants of the initial- 
value problem of the NLS equation [4]. The same result applies to periodic 
boundary-value problems. For finite line problems subject to homogeneous 
Dirichlet or Neumann boundary conditions at both boundaries, Eq. (5) 
indicates that the total mass is conserved, whereas the total linear momentum 
(cf. Eq. (6)) is not, in general, conserved. 

Similar conservation equations to those derived above can be obtained in 
nth-dimensional space (n=2,3) where the NLS equation can be written in 
dimensionless form as 

im = -V 2 u - \u\ 2 u, xGP, £>0. (8) 

In nth-dimensional space, the mass, linear momentum and energy densities 
are defined as 

m = u*u, M = i (uVu* - m*Vm) , E = Vu* ■ Vu - i |n| 4 , (9) 

where the mass and energy densities are scalars and the momentum is a vec- 
tor. Using Eq. ([8]), one can easily obtain the following conservation equations 

^ + V-M = 0, (10) 

= -V (4Vw*-V« - V 2 m - m 2 ) = -V (m 2 + AE - V 2 m) , (11) 
at v ' v ' 



-q1^ u * + -Qf Vu J = iV ' (Vm*V 2 m - VuVV) + V-(mM), 



(12) 

which indicate that the total mass, momentum and energy are conserved, 
i.e., they are invariant, for the initial- value problem of the nth-dimensional 
NLS equation if \u\ — > as |x| — > oo. Note that Eqs. (10)-(12) reduce to 
Eqs. (5)-(7) for one-dimensional problems. 

From the conservation of linear momentum (cf. Eq. (11)), i.e., from 
Newton's second law of classical mechanics, one can easily determine the 
classical potential and force densities as 

V d = m 2 + AE- V 2 m, f d = - W d , (13) 

where V c i and f c / denote the classical potential and classical force densities, 
respectively. Furthermore, one can easily deduce that the classical potential 
and force densities (cf. Eq. (13)) are related to the quantum energy and force 
densities (cf. Eqs (2) and (3) in nth-dimensional space) as 

e q = \{V cl - V 2 m-3m 2 ), (14) 

f d = 4f, - V(V 2 m + 3m 2 ) = f n + 2f d - 2V (Vu*-Vu) . (15) 

It is interesting to note that the above conservation equations for the 
mass, linear momentum and energy densities have been derived from the NLS 
equation, and have identical form to those of classical fluid dynamics [5] and 
electrodynamics [6,7]. Note that, in fluid mechanics, the mass conservation 
equation is also referred to as Lavoisier's law or continuity equation, while, 
in electrodynamics, the conservation of charge is referred to as continuity 
equation. 
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It is also remarkable that the mass and linear momentum equations de- 
rived in this paper have the same form as those of a lossless, nonlinear trans- 
mission line [8], i.e., 

dv(x,t) T dl{x,t) dl(x,t) dv{x,t) 

-^x- = - L ^^' -^ x - = - c{v) ^r~ (16) 

where v and / are the voltage and current, respectively, and L and C(v) 
denote the inductance and capacitance of the nonlinear transmission line per 
unit length, respectively. In this case, m, M and V c \ may be identified with 
the capacitance charge, Q, I and v/L, respectively. Furthermore, nonlinear 
LC ladder networks [9] may result in evolution equations of the Korteweg- 
de Vries (KdV) type which have soliton solutions, and asymptotic methods 
permit to deduce the NLS equation from the KdV equation [10]. 

3 Force Densities on Solitons in Infinite, Non- 
linear, Planar Waveguides 

In this section, the classical and quantum energy and force densities acting on 
a soliton propagating in an infinite, weakly nonlinear, planar waveguide are 
determined analytically. The analytical solution to the initial-value problem 
of Eq. Q can be written as 

u(x, t) = A sech£ exp ir) (17) 

where 

A 1 T / _ r 2 \ 1 

(18) 



t A ( A 1 

£= -j={x-X Q -Ct), T]=- 



<•(.!• - .r„) ■:■ {A 2 - C -)l O u 
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where A, c, x and O are the soliton's amplitude, speed, initial position and 
intial phase, respectively. 

Using this solution, one can easily obtain the following quantum densities 

e q =(j~^-)A^ech^ (19) 



f n = -2V2A 5 sech 4 ^ tanh f (20) 

Rprh 2 
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A3 r -1 

/ d = — sech 2 £ tanh f (c 2 - 2A 2 ) + Sv^A 2 sech 2 £ (21) 



f q = fn + fd = ~^-A 3 (c 2 - 2A 2 ) sech 2 £ tanh £, (22) 
and the following classical potential and force densities 

V d = A 2 c 2 sech 2 £ (23) 

f d = V2A 3 c 2 sech 2 £ tanh f . (24) 

From Eqs. (19) and (23) and Eqs. (22) and (24), the following expression 
can be readily obtained 

V cl -f cl ~4 2c 2 ' 

which indicates that, for the initial-value problem of the NLS equation, the 
quantum and classical potential and force densities are proportional to their 
classical counterparts. In fact, the magnitude of the quantum force density 
is 25 per cent of that of the classical one for A—c—1 as shown in Figure 
1. Note that the difference in sign between the classical and quantum force 
densities for A—c—1 is due to the definitions of the quantum energy and 
classical potential (cf. Eqs. (2) and (13)). 
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4 Force Densities on Solitons in Finite, Non- 
linear, Planar Waveguides 

In finite, nonlinear, planar waveguides, Eq. (1) was solved numerically by 
means of a second-order accurate, in both space and time, Crank-Nicolson 
method subject to 

— ( x ,t) = -j— (x + 2kH,t), Vn > 0, k € Z, x E V, t > 0, (26) 

for periodic waveguides, and 

u(-H,t)+'yu x (-H,t) = 0, u(H,t)+'yu x (H,t) = 0, t > 0, (27) 

where V=[—H, H], H = 50, and the values 7=0, oo and 7 ^ correspond to 
homogeneous Dirichlet, Neumann and Robin boundary conditions, respec- 
tively, at both boundaries of the waveguide. 

The time step and the grid size used in the calculations shown in the next 
section were 0.01 and 0.25, respectively, </)o=Xq=0 and A—c—1. The initial 
condition corresponded to the exact solution of the initial-value problem of 
the NLS equation, i.e., Eq. (17), translated in such a manner that there 
were not mathematical incompatibilities between the initial and boundary 
conditions. 

5 Presentation of Results 

Figures 1-3 show the classical and quantum force densities acting on solitons 
which propagate in weakly nonlinear, planar waveguides. The left side of 
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Figure 1 shows the force densities in an infinite waveguide, i.e., those cor- 
responding to the initial-value problem of the NLS equation, and indicates 
that the magnitude of the classical force density is four times larger than the 
quantum one in agreement with the analytical results presented in Section 3. 
The left side of Figure 1 also shows the difference in sign between the classical 
and quantum force density which is due to the definitions of quantum energy 
and classical potential employed in this paper, and the S'-shape of both force 
densities. 

The results presented in the right side of Figure 1 correspond to a finite 
waveguide subject to periodic boundary conditions. Except for truncation 
errors, the classical and quantum force densities presented in the right side of 
Figure 1 coincide with those of the left side of the same figure. Furthermore, 
the linear momentum does not change as the soliton interacts with the right 
boundary in agreement with the analysis of Section 2. 

For a finite waveguide subject to homogeneous Dirichlet boundary condi- 
tions at both boundaries, the left side of Figure 2 indicates that the classical 
force densities prior to and after the collision of the soliton with the right 
boundary are identical, and that, just during the soliton's interaction with 
the boundary, the classical force density undergoes large changes at the right 
boundary. The maximum and minimum values of the classical force density 
are 0.5637 and -2.2249, respectively, whereas those of the quantum force den- 
sity are 0.8059 and -0.7825, respectively. Note that, prior to and after the 
interaction of the soliton with the right boundary, the largest values of the 
classical and quantum force densities are 0.5443 and 0.1361, respectively, as 
shown in the upper left corner of Figure 1. 
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The left side of Figure 2 shows the classical and quantum force densities 
on solitons propagating in finite, planar waveguides subject to homogeneous 
Neumann boundary conditions at both boundaries. Due to the large changes 
in both the classical and the quantum force densities during the collision of 
the soliton with the right boundary, the S'-shape of these force densities prior 
to and after the collision cannot be observed in the right side of Figure 2. 
The maximum and minimum values of the classical force density are 1.1283 
and -133.2461, respectively, whereas those of the quantum force density are 
6.9607 and -44.4007, respectively. 

The results presented in Figure 2 clearly indicate that the magnitude of 
the classical force density is, in general, much larger and undergoes greater 
changes than the quantum one during the interaction of the soliton with the 
right boundary. 

Since the one-dimensional NLS equation subject to homogeneous Dirich- 
let or Neumann boundary conditions is invariant under mirror reflections in 
x, the interaction of the soliton with the left boundary is expected to be 
identical to that with the right one. For this reason, only the interaction of 
the soliton with the right boundary was shown in Figure 2. 

The results presented in Figure 3 correspond to a finite, planar waveguide 
subject to homogeneous Robin boundary conditions at both boundaries and 
7=1 (cf. Eq. (27)). Since the one-dimensional NLS equation is invariant 
under reflections in x, while the homogeneous Robin boundary conditions 
are not, the interaction of the soliton with the right boundary is expected to 
be different from that with the left one as illustrated in the left and right, 
respectively, sides of Figure 3. This figure clearly indicates that, due to 
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the large increase that the classical force density undergoes as the soliton 
interacts with the right boundary, the S-shape of the classical force density 
can be barely seen. The maximum and minimum values of the classical 
force density during the interaction of the soliton with the right boundary 
are 21.7091 and -1.3303, respectively, whereas those of the quantum force 
density are 3.5286 and -0.7605, respectively. 

Figure 3 also indicates that the classical force density undergoes larger 
variations than the quantum one and larger variations as the soliton inter- 
acts with the right boundary than when it interacts with the left one. The 
maximum and minimum values of the classical force density during the inter- 
action of the soliton with the left boundary are 0.8590 and -6.1238, respec- 
tively, whereas those of the quantum force density are 0.5064 and -2.3679, 
respectively. 

6 Conclusions 

The nonlinear Schrodinger equation has been used to determine the conser- 
vation equations of mass, linear momentum and energy densities of solitons 
propagating in infinite or finite, weakly nonlinear, planar waveguides. It has 
been shown that these conservation equations have the same form as those 
of classical fluid dynamics and electrodynamics. It has also been shown 
that the classical force density is larger than the quantum one and that 
it undergoes large changes as the soliton interacts with boundaries subject 
to homogeneous boundary conditions. The largest changes in the classiccal 
force densities occur in planar waveguides subject to homogeneous Neumann 
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boundary conditions. 
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Figure 1 Nonlinear (/„), diffraction [Re(/ rf )] ( and total (/ tota ,) force 
densities and momentum density (p) for periodic boundary conditions 
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Figure 2 Nonlinear (/„), diffraction [Re(/ d )], and total (/, ota i) force 
densities and momentum density (p) for Dirichlet boundary condi- 
tions 
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Figure 3 Nonlinear (/„), diffraction [Re(/ d )], and total (/, ota i) force 
densities and momentum density (p) for Neumann boundary con- 
ditions 
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Figure 4 Nonlinear (/„), diffraction [Re(/J], and total (/total) force 
densities and momentum density (p) for Robin boundary conditions 
and first collision with the right boundary 




Figure 5 Nonlinear (/„), diffraction [Re(/ rf )], and total (/ tola i) force 
densities and momentum density (p) for Robin boundary conditions 
and first collision with the left boundary 



